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ABSTRACT 

Previous studies of the interior structure of transiting exoplanets have shown that the heavy element 
content of gas giants increases with host star metallicity. Since metal-poor planets are less dense and 
have larger radii than metal-rich planets of the same mass, one might expect that metal-poor stars 
host a higher proportion of gas giants with large radii than metal-rich stars. Here I present evidence 
for a negative correlation at the 2.3er level between eclipse depth and stellar metallicity in the Kepler 
gas giant candidates. Based on Kendall's r statistics, the probability that eclipse depth depends 
on star metallicity is 0.981. The correlation is consistent with planets orbiting low-metallicity stars 
being, on average, larger in comparison with their host stars than planets orbiting metal-rich stars. 
Furthermore, since metal-rich stars have smaller radii than metal-poor stars of the same mass and 
age, a uniform population of planets should show a rise in median eclipse depth with [M/H]. The 
fact that I find the opposite trend indicates that substantial changes in gas giant interior structure 
must accompany increasing [M/H]. I investigate whether the known scarcity of giant planets orbiting 
low-mass stars could masquerade as an eclipse depth-metallicity correlation, given the degeneracy 
between metallicity and temperature for cool stars in the Kepler Input Catalog. While the eclise 
depth-metallicity correlation is not yet on firm statistical footing and will require spectroscopic [Fe/H] 
measurements for validation, it is an intriguing window into how the interior structure of planets and 
even the planet formation mechanism may be changing with Galactic chemical evolution. 
Subject headings: Stars: planetary systems — Planets and satellites: formation — Planets and satel- 
lites: fundamental parameters — Planets and satellites: composition 



1. INTRODUCTION 

Since 1958, astronomers have known that Jupiter and 
Saturn are smaller than they wo uld be if they we re made 
of pure hydrogen and helium ([Demarcusl 11958ft . Even 
though planetary radii depend on complex equations of 
state, the simple principle that metal-rich planets are 
more dense than their metal-poor counterparts of the 
same mass holds for a wide r ange of planetary radii (e.g. 
iFortnev & Nettelmanr] 120111) . Since the heavy-element 
content of t ransiting gas giants increases with host sta r 
metallicity (|Guillot et al.l 120061: iMilTer fc Fortnevl 120111 ). 
one may make a simple prediction: the gas giants orbit- 
ing metal-poor stars should have higher radii, on aver- 
age, than the gas giants orbiting metal-rich stars. Just 
as the planet detection rate depends on star metallic- 
ity dGonzalezlll998t [Santos et al.|[200l iFischer fc Valentil 
12005ft . the typical giant planet interior structure also 
should change according to the metal content of the host 
star. 

Before the advent of the Kepler mission, there was no 
large, uniform sample of gas giants available to search 
for a relationship between stellar metallicity and gas 
giant size. The surveys of transiting gas giants were 
heterogeneous and mostly insensitive to planets smaller 
than ~ 0.7i?j up (0.7 Jupiter radii). One notable excep- 
tion was GJ 436 b at 0. 365-Rj UD , wh i ch orb its an M2.5 
star of radius 0.464i? Q (jGillon et alj l2007bf) . The Ke- 
pler spacecraft, launched in March of 2009, is monitoring 
about 150,000 mostly Solar-type stars, primarily search- 
ing for transits of Earthlike planets in the habitable zone 
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(jBatalha et a,l.|[20Toh . Yet the Kepler mission is also ex- 
traordinarily useful for studies of gas giants, for the first 
time providing a sample of candidate giant planet hosts 
observed under uniform conditions and with the sensitiv- 
ity to detect all of the Neptune-like objects of ~ 0.3i?j up 
and below. 

Here, then, is the opportunity to investigate the de- 
pendence of planetary interior structure on host star 
metallicity. Of course, one expects to see gas giants 
with a range of sizes orbiting all types of Population 
1 stars — after all, Jupiter and Neptune orbit the same 
star but have radii that differ by almost a factor of 
three. However, the ensemble of planet radii should 
show some dependence on stellar metallicity if the heavy 
element contents of planets an d their paren t star s are 
truly connected, as sugg ested bv lGuillot et al.l (|2006l ) and 
Millcr "fc Fortnevl (|201lD . Likewise, since the Kepler mis- 
sion is surveying primarily similar, Sunlike stars, one also 
expects a negative correlation between gas giant eclipse 
depth and star metallicity. In this article I present the 
first evidence that such a trend exists — though at this 
time the eclipse depth-metallicity trend is on tenuous sta- 
tistical footing and will require follow-up spectroscopic 
observations for verification. 

This article is organized as follows: In <J2]I outline my 
gas giant selection criteria. In SJ3]I present statistical sup- 
port for a negative correlation between eclipse depth and 
star metallicity. I discuss possible selection biases that 
could masquerade as an eclipse depth-metallicity trend in 
||4] Finally, I examine physical mechanisms besides the 
simple density effect that could lead to an eclipse depth- 
metallicity correlation and speculate about the implica- 
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tions for planet formation theory in [JSJ 
2. SAMPLE SELECTION 

In February 2011, the Kepler team released data for 
1235 planet candidates ([Borucki et all 1201 lal) . While 
the planet/star radius ratios R p /R* were directly mea- 
sured from eclipse depths, the values of i?» were based 
on log(<7) measurements from the Kepler Input Cata- 
log (KIC: iBrown et al.|[20Tl . The KIC stellar param- 
eters were computed using photometry in the Sloan u, 
g, r, i, z filters, the intermedia te-bandwidth D51 fi l- 
ter, and the 2 MASS JHK filters (jSkrutskie et al.ll200fih . 
IBrown et al.l (|2011[) quote a 0.4-dex error bar on log(g), 
which translates into a +58%/-37% radius error for a 
planet transiting a Sunlike star. Hotter, early- to mid- 
F stars have higher log(g) errors, while subgiants tend 
to ha ve systematica l ly hig h log(g) estimates. Further- 
more, iChaplin et al.l ([20101 ) performed a detailed astero- 
seismic follow-up study of three Sunlike stars and found 
that the KIC T e g estim ates were systematically low. 
iPinsonneault et all (|2012f ) also found that temperatures 
derived from KIC griz photometry were systematically 
lower than tempera tures computed by appl ying the in- 
frared flux method (iCasaerande et~atl 120101 ) to 2MASS 
photometry (jSkrutskie et al.ll2006f) . The problems with 
log(g) for hot stars may therefore apply to stars that have 
apparently Sunlike temperatures according to the KIC. 
Quoted planet candidate radii from the Kepler data re- 
lease are subject to significant uncertainties. 

However, the Kepler planet candidate radii can still 
be used to identify a set of likely gaseous, giant planet 
candidates. Here I define giant planets as objects with 
bulk density p < 1.5 g cm -3 , for which ~ 50% or 
more of the planetary volume is due to the H/He at- 
mosp here. According to the planetary interior mod- 
els of Helled ct al. (20ll|), this volume-based definition 
of giant planet would include Uran us and Neptune, 
which are about 65% H/He by v olume. iRogers fe Seagerl 
(|2010f) sh ow that HAT-P-llb (iBakos et all 12010ft and 
GJ 436b ([Gillon et al.l l2007at iTorresI l2007n must also 
have gas mass and volume fractions similar to Uranus 
and Neptune. HAT-P-llb, GJ 436b, Uranus and Nep- 
tune have radii of 4.7i?®, 4.2i?®, 4.1i?®, and 4.0i?®, 
respectively, suggesting that gas giants in general have 
radii R n > 4J? m . To be conservative, I follow 
ISchlaufman fc Laughlinl ()2011l ) and set the minimum 
planet candidate radius for this study at R p = 5.0i?®. 
For a 5i?® candidate with a +58% measurement er- 
ror on planet size, the true size would be 3.2i?®, still 
within the size range wh ere a gaseous atmosphere is likely 
(|Rogers fc Seagerl [20101) . 

The basic sample consists (1) of candidate giant plan- 
ets with radius R p > 5.0i?® that (2) also have metallicity 
measurements from the KIC. Of the candidates that sat- 
isfy criteria (1) and (2), 28 have estimated sizes larger 
than the "styrofoam planet" Kepler-7b {R p — 18. 1-R®; 
lLatham et al.ll20i0ft , the lowest-density transiting planet 
yet known. Given that theoretical mass-radius relations 
predict that Jupiter-mass (Mj up ) planets, 20Mj up brown 
dwarfs a nd 100M,T„p stars hay e radii that vary by less 
than 5% (|Chabrier et al.l [20091) . it is unfortunately very 
easy for eclipsing low-mass companion stars to masquer- 
ade as giant planets. I therefore carefully scrutinized 
the light curves of the largest planet candidates for signs 



of (a) planet radius mismeasurement due to V-shaped 
eclipses, (b) odd/even eclipse depth differences that in- 
dicate eclipsing binaries (either in the target itself or a 
background object), and (c) pulsations that suggest the 
8 Scuti-like variability of early- to mid-F stars, for which 
the \og{g) uncertainties are high. 

Among the 23 candidates with metallicities from KIC 
and R p > 18. li?®, I find 7 objects for which the re- 
ported ratio R P I R* was mis-fit by the Kepler data, analy- 
sis pipeline. All such light curves show V-shaped eclipses 
rather than flat-bottomed eclipses, which drives the fit- 
ting algorithm toward artificially low impact parameters 
and high values of R p / R„ (B. Cochran, private commu- 
nication). Two objects, KIC 5649956 and KIC 1432214, 
show the odd/even differences in eclipse depth that indi- 
cate an eclipsing binary system. KIC 5649956 pulsates 
at the 3-4% level, making the odd/even effect difficult 
to identify. KIC 6470149 and KIC 7449844, which do 
not have odd/even eclipse patterns, have been shown to 
host eclipsing M dwarfs by follow-up RV observations 
taken since the release of the Kepler planet candidate 
list. Meanwhile, follow-up spectra of KIC 5356593 have 
raised its log g value from the 3.8 dex listed in the KIC to 
4.5 dex, correspondingly shrinking the computed planet 
radius from 35. 5i?® to 13-14i?® (B. Cochran, private 
communication) . 

Furthering the list of suspicious objects, six other stars 
have no obvious mismeasurement of eclipse depth, but 
their eclipses are V-shaped, which is more consistent with 
a companion star than a planet. One star, KIC 10616571, 
pulsates at the 3% level on three frequencies, while one 
or two small peaks occur at the minimum of each eclipse, 
suggesting a hierarchical triple system. Finally, five stars 
that do not have V-shaped eclipses or odd/even effects 
pulsate with similar amplitude to the transit signal. One 
such pulsator, KIC 11818800, will have its planet can- 
didate radius revised from 36i?® to 9.77i?® in the up- 
coming re- release of the Kepler planet candidate list (B. 
Cochran, private communication). The preponderance 
of problems among the largest planet candidates suggest 
that one should exercise caution when selecting a sample 
of inflated giant planets. 

In a planet candidate list ordered according to R p , 
I found that pulsations, odd/even effects and mismea- 
sured and/or V-shaped eclipses accounted for fully 100% 
of the objects with R p > 22i?®. Below 20i?®, light 
curves that have no obvious problems begin to domi- 
nate. The confirmed planet Kepler-12b has R p = 18. 9i?® 
(|Fortnev et al.ll20TU) and a similar density to Kepler- 7b, 
validating the idea that such inflated planets can exist. 
I therefore set an upper limit of 20i?® to the planet can- 
didates in this study. The basic sample consists of 213 
giant planet candidates with 5.0-R® < R p < 20i?® whose 
host stars also have metallicities from the Kepler Input 
Catalog (KIC). All planet candidates used in this study 
are liste d in Table 1. 

While IMorton fc Johnson! (120111) argue for a false pos- 
itive rate of 5%-10% amon g the Kepler planet candi- 
dates, iBOTUckTeTliI] <[20TT3) are more cautious and quote 
> 80% reliability for rank 2 KOIs (Kepler Objects of In- 
terest) and > 60% reliability for rank 3 and 4 KOIs. De- 
spite the scrutiny of planet candidates, there is therefore 
still the possibility that the sample is contaminated by 
as much as 40% eclipsing binaries, background eclipsing 
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binaries and hierarchical triples. I discuss the possible 
consequences of such contamination in i|4j 

3. STATISTICAL EVIDENCE FOR DEEP ECLIPSES OF 
LOW-METALLICITY STARS 

With a sample of gas giant candidates in hand, I now 
present the first statistical evidence that Galactic chem- 
ical evolution affects the internal structure of giant plan- 
ets as well as the frequency with which they form. Figure 
Q] shows R p /R* as a function of host star [M/H] for the 
stars in my sample. Planet candidates are color-coded 
according to the irradiation they receive from their host 
stars, which I discuss further in fj5] Although I used 
the Kepler estimates of R p to select my sample of gas 
giant candidates (see for selection criteria) , for the 
purposes of statistical analysis I have cast the trend in 
terms of R p /R* rather than R p because planet radii 
tend to change substantially between inclusion in the 
Kepler planet candidate list and eventual confirmation 
and publication. R p /R*, by contrast, is a simple observ- 
able whose accuracy I have verified by examining the 
light curves of the largest planet candidates. The ques- 
tion I am asking is — are Kepler candidate giant plan- 
ets orbiting low-metallicity stars larger in comparison 
to their host stars than candidate planets orbiting high- 
metallicity stars? 

Figure [1] shows eclipse depth as a function of [M/H] 
for the candidate giant planet hosts in the selected sam- 
ple. Scanning from left to right, notice the knot of small 
planet candidates with 0.05 < R p /R* < 0.13 that ap- 
pears at [M/H] > —0.2 without a commensurate increase 
in large planet candidates with R p /R* > 0.13. To quan- 
tify the trend in decreasing eclipse de pth with metallic- 
ity, I have computed a running median (Hard le fc Steigeil 
119951) with a 21-point window (solid black) and a local 
polynomial regression (Clevelan d et al.lll992[ ) with span 
0.75 and degree 1 (dash-dot black). Notice the decreasing 
nature of both manifestations of the trend in the range 
-0.1 < [M/H] < 0.3. While planets of all radii can ap- 
pear around stars of all metallicities, there is a hint that 
eclipses may be deepest, on average, for low-metallicity 
stars. This possible trend is as predicted — planets orbit- 
ing low-metallicity hosts should have the lowest solid/gas 
ratios, making them the least dense and therefore biggest 
of planets. 

Given that the running median is noisy and the slope 
of the local polynomial regression is small, is the trend 
statistically significant? The Kendall rank coefficient for 
mapping [M/H] onto R p /R* is r = —0.108. The associ- 
ated probability that the null hypothesis — that R p / i?* is 
unrelated to star metallicity — is correct is p = 0.019, a 
significance level at which the R statistical suite recom- 
mends adoption of the alternative hypothe sis. F o r mor e 
information on Kendall's r coefficient, see lAbdl (f2007h . 
Using the large-sample Gaussian approximation to the 
variance v of the r sampling distribution, 



2 (2n + 5) 
9n (n — 1) ' 



(1) 



where n = 213 is the number of stars in the sample 
(|Abdil [20071 ). I find that the standard deviation of the 
t sampling distribution is 0.046. The computed value of 
t = —0.108 differs from the r = expected in the case 
of the null hypothesis by — 2.3er. I therefore consider 



the eclipse depth-metallicity trend interesting and sug- 
gestive, but not conclusive — especially since so little is 
known about the true distribution of giant planet eclipse 
depths that the Gaussian approximation to the t distri- 
bution may not be appropriate. My main goal in writing 
this article is to motivate spectroscopic follow-up of the 
Kepler giant planet candidates. More precise metallici- 
ties and star radii will place our understanding of how 
planet structure evolves with star metallicity on much 
firmer statistical footing. 

Besides the robust measures of center polotted in Fig- 
ured] and the non-parametric Kendall's r coefficient dis- 
cussed above, there is another way to visualize the eclipse 
depth-metallicity trend. In Figure [2] I have divided the 
data into four broad metallicity bins and plotted the 
fraction of planet candidates for which R p /R* > 0.13 
in each bin. I chose R p /R* = 0.13 as the cutoff be- 
cause it is a roughly typical eclipse depth of well-studied 
inflated planets such as HP 209458 b (R p /R* = 0.12; 
iCharbonneau et alj [2000h and HAT-P-19 b (R p /R* = 
0.14; lHartman et al.ll20"Tll) . What Figure [2] also hints at, 
then, is that stars in the lowest-metallicity bin tend to 
host a higher fraction of planets that are large in com- 
parison to their host stars than the stars in the highcr- 
metallicity bins. Notice, however, that the Poisson error 
bars on the fraction of high eclipse depths in the lowest- 
metallicity bin overlap with the error bars on two of the 
other three bins. Furthermore, note that in binning the 
data and setting a threshold value of R p /R* — 0.13, 
I have made the threshold and the bins parameters 
that affect any statistical inferences made from Figure 
[2] (though slightly different cutoff values yield a similar 
trend). Figure [2] is intended more as a visualization of 
the data than as a statistical tool. 

There is one more potential pitfall left to examine, 
the accuracy of the KIC [M/H] estimates. Already the 
planet size-metallicity trend is on marginal footing — 
could [M/H] inaccuracies wipe it out entirely? Figure 
[3] shows KIC [M/H] as a function of spectroscopically 
determined [Fe/H0 for confirmed, published Kepler dis- 
coveries. From Figure [31 one can see that the KIC metal- 
licities are systematically lower than the spectroscopic 
values. The best-fit line relating KIC [M/H] to spectro- 
scopic [Fe/H] is 

[M/H] K ic = 0.88[Fe/H] spcc - 0.14. (2) 

This result agrees with the analysis of iBrown et al.1 
(pOll . whose comparison of KIC metallicities with 
[Fc/H] derived from Keck/HIRES spectra uncovered 
a systematic 0.17-dex underestimate in the KIC val- 
ues. The Spearman rank correlation coefficient relat- 
ing [Fe/H] spcc and [M/H]kic is 0.76, which has a two- 
sided signifi cance of its deviatio n from zero of2.7xl0~ 6 . 
In contrast. IBrown et alJ (|2011l ) found a Spearman rank 
correlation coefficient of only 0.42, with a significance 
of 0.02, for a subset of stars in the KIC observed with 
Keck/HIRES. My analysis indicates that, despite the dif- 
ferences in absolute scale between spectroscopic and KIC 
metallicities, the two are nevertheless closely related. 
Note, however, the two outliers, Kepler-10 and Kepler- 
19. Kepler-10 is a Solar-radius, Solar-temperature star 

1 Host star metallicities for Kepler discoveries were collected 
from the Extrasolar Planets Encyclopaedia, http://exoplanet.eu 
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with G spectral type (|Batalha et al.ll2011[) , while Kepler- 
19 is a slightly smaller star (i?* = 0.85Rr^. [Ballard et al.l 
1201 ID with T eff = 5541 ± 60 K and a possibly later spec- 
tral type. Neither star falls into a known problem cat- 
egory for the KIC — F stars, subgiants and stars with 
T c ff < 4200 K — so there is no clear reason why their 
KIC metallicities are so low. Published planet candi- 
dates used in the [M/H] accuracy analysis are listed in 
Table [2] 

Knowing that KIC metallicities do track spectroscopic 
metallicities, albeit with some scatter, the relevant ques- 
tion for the purposes of this study is whether or not one 
can place the Kepler candidate planet hosts in the cor- 
rect order on some metallicity scale. Keeping in mind 
the warnings about how binning introduces extra param- 
eters into the statistical analysis whose effects may not 
be well understood, I still wish to know how often mea- 
surement error will cause a candidate planet host to cross 
over into a neighboring bin in Figure [5J The scatter of 
the KIC [M/H] values around the best- fit line in equa- 
tion [2] is C[m/h] = 0.22 dex, indicated by dash-dotted 
lines in the left-hand panel of Figure [3] (The right-hand 
panel of Figure [3] shows KIC [M/H] as a function of T c g , 
a correlation I will discuss further in SJU) I modeled the 
probability distribution of each star's true metallicity us- 
ing a Gaussian distribution with a = 0.22 dex centered 
on that star's KIC [M/H] value. These Gaussian dis- 
tributions of probable metallicity formed the basis of a 
Monte Carlo simulation, in which I randomly sampled 
each star's [M/H] distribution 20,000 times, divided the 
random samples into the same [M/H] bins as in Figure 

and counted the fraction of planet candidates in each 
bin with FL/R* > 0.13. 

Figure [5] shows the results of the Monte Carlo sim- 
ulation. Here the story from the error bars is slightly 
more promising. Candidate planet hosts with —0.85 < 
[M/H] < —0.45 dex show significantly more deep eclipses 
for which R p /R* > 0.13 than stars with [Fe/H] > 
—0.15 dex. Once again, though, the error bars barely 
miss overlapping, so that the statistical significance of 
the trend is about 2<r — suggestive, but not definitive. 
Confirming that gas giants orbiting low-metallicity stars 
really tend to have larger radii than gas giants orbit- 
ing high-metallicity stars will require extensive follow-up 
spectroscopy of the Kepler candidate planet hosts. The 
California-Kepler survey is already underway and will 
significantly impr ove stellar log(q), [F e/H] and planet 
radius estimates (|Howard et al J 120121 ) . With spectro- 
scopic observations of a large fraction of the giant planet 
candidates, it should even be possible to cast the size- 
metallicity trend in terms of R p rather than R p /R*. 

4. STATISTICAL BIASES 

In this section I discuss possible statistical biases that 
could masquerade as a size-metallicity trend. The first 
possible bias I wish to examine is whether or not the Ke- 
pler mission is equally sensitive to planets with minimum 
radius 5i?® in each metallicity bin. Given its ability to 
detect Earthlike planets transiting Sunlike stars, Kepler 
has at least 25 times the sensitivity required to detect 
planets with R p — 5i?® orbiting Solar-type stars, which 
produce eclipse depths of two parts per thousand. Plan- 
ets with R p — 5Rq only become undetectable when the 
star radius approaches 5Rq. Since the largest star radius 



in the selected sample is 4.82i? Q and only five stars have 
i?* > 3i?0 , there should be no systematic trend that pre- 
vents detection of the smaller, 5i?0 giant planets at low 
metallicity. 

Another possible bias arises from the interaction of two 
effec ts: the scarcity of g as giants orbiting low-mass stars 
(e.g. lBonfils et al.ll2011l ) and the red colors of late K- and 
M-type stars. Among stars of a constant J — H color, 
those that are most metal-rich tend to have the highest 
g — r color because of the wealt h of iron absorption lines 
in the blue part of the spectrum (Schlauf man fc Lau ghlin 
1201 If ) . Red optical color is therefore an important metal- 
licity indicator for stars in the KIC, but its usefulness 
breaks down for late-type st ars that are natural ly red, 
whatever their metallicity. iBrown et al.l (|2011| ) show 
that the coolest stars, those with T c g < 4200 K, tend 
to be artificially classified with super-Solar metallicity. 
Above T c g = 4200 K, the metallicity distribution is a 
well-behaved Gaussian that is not a function of tempera- 
ture. One might worry that the low proportion of inflated 
giant planets orbiting the most metal-rich stars could 
be an artifact of the known scarcity of gas giants or- 
biting low- mass stars (|Endl et all 120061 : iCumming et al.l 
l2008tlBonfils et alJl2011[ ). Perhaps low -mass stars whose 
metallicities are artificially high account for much of the 
apparently metal-rich part of the selected sample. For- 
tunately, only four stars in the sample have T c ff < 4200, 
so the most problematic KIC metallicities have largely 
been avoided. 

Altho ugh the ana l ysis o f KIC atmospheric parame- 
ters by IBrown et al.1 (|2011h does not show any obvious 
correlation between KIC [M/H] and T c g for stars with 
T c ff > 4200 (see their Figure 11), I wish to independently 
verify the lack of a temperature-metallicity relationship. 
The right-hand panel of Figure [3] shows [M/H] as a func- 
tion of T e ff for stars in the selected sample. The solid 
line indicates the mean [M/H] value and the dash-dot 
line gives the best linear fit. The Spearman rank corre- 
lation coefficient is -0.19, with a significance of 0.0049 — 
three orders of magnitude higher than the significance 
of the [M/H] Kic- [Fe/H] spcc relationship (a high signifi- 
cance value indicates a weak correlation). While the cor- 
relation analysis shows some evidence for a metallicity- 
temperature relationship, KIC metallicity estimates ap- 
pear to be much more related to true star metallicities 
than temperatures. It is likely that the eclipse depth- 
[M/H] correlation heralds a change in planet structure 
with host star [M/H] rather than with host star mass 
and temperature, but I am unwilling to completely rule 
out the possibility that dearth of gas giants orbiting low- 
mass stars is affecting the results. 

Although the Kepler planet radii used to select the 
sample show no correlation with log(g) (Figure |4j left 
panel), a possible degeneracy between [M/H] and log(p) 
in the KIC could affect the results of this study. Red 
optical color is an indicator of both high log(g) and high 
metallicity, and the only intermediate-bandwidth filter 
used for KIC photometry — the T>51 filter — is sensitive 
to both log(<?) and [M/H]. Fortunately, the correlation 
between log(g) and [M/H] is much weaker than the cor- 
relation between T e g and [M/H]. For [M/H] as a func- 
tion of log(g) (Figure |H right panel), I find a Spearman 
rank correlation coefficient of -0.13, with a significance 
of 0.05. Furthermore, the subgiant problem discussed 
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earlier in this section should apply to a small fraction 
of the stars in the sample: only seven stars of 213, or 
3%, have log(g) < 0.7. As with temperature, I am un- 
willing to definitively state that the weak degeneracy 
between KIC log(g) and [M/H] plays no part in creat- 
ing the eclipse depth-metallicity correlation. However, I 
think it likely that declining eclipse depths with metal- 
licity are related to giant planet interior structure rather 
than log(<?)-[M/H] degeneracies. 

The final possible source of bias comes from the possi- 
ble 5%-40% contamination of the sample with eclipsing 
binaries, background eclipsing binaries, and hierarchical 
triples. This bias arises because stellar multiplicity rates 
depend on star mass. From O.IMq — 40M Q , the bi- 
nary fraction is a n increasing function of star mass (e.g. 
iClark et al.ll2012t ). with about 20% of O.1M stars in bi- 
naries and 40% of IMq stars in binaries. The increasing 
binary fraction with Af* may mean that the planet candi- 
dates orbiting stars with the largest radii have the highest 
contamination rate. Since the tendency in the KIC is for 
low-mass, low-temperature stars to have artificially high 
metallicity values, it is possible for background eclipsing 
binaries to produce an eclipse depth-metallicity trend. 
However, the binary fraction is only a source of bias if 
[M/H] is related to star temperature in the candidate gi- 
ant planet sample. Given that my analysis indicates that 
KIC metallicity is more related to spectroscopic metallic- 
ity than any other stellar parameter, I do not think the 
multiple-star contamination is the source of the possible 
eclipse depth-metallicity trend — but, again, more follow- 
up is necessary to rule out multiple stars as a source of 
bias. 

5. PHYSICAL EXPLANATIONS FOR THE ECLIPSE 
DEPTH-METALLICITY CORRELATION 

Assuming that the tentative correlation between planet 
candidate eclipse depth and metallicity (1) is real, and 
(2) is not a manifestation of the dependence of planet 
occurrence rate on stellar mass — both of which are as- 
sumptions that require further investigation — what is 
the physical reason for the correlation? First, it is 
important to point out that for stars of a given mass 
and age, s tar radius decrease s with increasing metal- 
licity (e.g. iMarigo et al.l |2008|) — so the fall in median 
Rp/R* with [M/H] means gas giant interior structure 
must change substantially to counteract the fact that 
eclipse depth should tend to rise with star metallicity 
for a uniform population of planets. To the extent that 
the KIC [M/H] values track spectroscopic [Fe/H], the 
eclipse depth-metallicity correlation reflects real changes 
in the population of giant planets as a function of star 
metallicity. 

One possibility is that irradiation from the cen- 
tral star depends on metallicity. Irradiation pro- 
vides the energy required to inflate hot Jupiters ei- 
ther through Ohmic dissipation (which depends on 
the star's ability to ionize the upper atmosphere; 
iBatygm et al.l 120111 ) or simply by forcing a shallow at- 
mospheric temp erature gradient that slows planetary 
contraction (e.g.lGuillot fc Showman 2002; Baraffe et al.l 
120031: iBurrows et al.ll2003l: IFortnev fe Nettelmannll2011[ ). 
At a given mass, a star of low metallicity will be slightly 
bluer and hotter than its counterpart at high metallic- 
ity, so metal-poor stars might tend to host larger-radius 



planets if the orbital distribution does not depend on 
metallicity. The points in Figure [1] are color-coded ac- 
cording to average stellar insolation. Although there is a 
visible tendency for least-irradiated planets to have the 
smallest eclipse depths (note the cluster of green points 
toward the bottom of the plot), there is no obvious ten- 
dency for the planets orbiting the lowest-metallicity hosts 
to be the most irradiated. I therefore tentatively rule 
out differences in irradiation levels as the reason for the 
eclipse depth-metallicity correlation. 

Other planetary inflation mecha nisms that have bee n 
proposed includ e tidal heating (Ijackson et ail [2008), 
thermal tides (lArras fc Socrates] 20101) and double- 
diffusive convection (jChabrier fc Baraffd f2007h . Of 
these alternative inflation mechanisms, double-diffusive 
convection — in which composition gradients create mul- 
tiple, semi-detached convective layers between which en- 
ergy transport is inefficient — is the most likely to be re- 
lated to star metallicity. In their mod els of the interior 
struct ures of uninflated hot Jupiters, iMiller fc Fortnevl 
(|2011[) found a positive correlation between stellar metal- 
licity and heavy-element mass in planetary envelopes 
(see their Figure 2). A metal- rich planetary envelope 
would likely have a steeper composition gradient than a 
metal-poor envelope. It is therefore possible that double- 
diffusive convection provides the required connection be- 
tween giant planet structure and stellar metallicity. 

The most straightforward explanation for the eclipse 
depth-metallicity trend, however, is that metal-rich 
planets of a given mass are denser than their metal- 
poor counterparts, leading to smaller radii (e.g. 
IFortnev fc Nettelmannl 1201 ID . Transit searches have 
for some time provided hints that heavily metal- 
enriched, dense plane ts tend to orbit metal- rich stars 
such as HP 1490 2 6 (ISato et al.l (2005^ and HAT-P-2 

and the con- 



(jBakos et al.l 120071 : iLeconte et aD 



2009) 



nection between heavy-element mass and host star 
metallicity — the crucial ingredient behind any eclipse 
depth- [M / H] correlation— has also been statistically es - 
tablished (|Guillot et al.l 120061 : IMiller fc Fortnevl 120 111 ). 
Given that Jupiter, Saturn, Uranus and Neptune all or- 
bit the same star but have radii ranging from A.ORq to 
11.2.R®, one expects to see gas giants with a range of radii 
orbiting stars of all metallicities yet probed (whether or 
not giant planets should form around Population II stars 
is an interesting question). However, if Jupiter, Saturn, 
Uranus and Neptune were all twice as metal-rich as they 
are, the average radius of the entire ensemble of Solar 
System planets would decline. I favor a simple density 
effect as the reason behind the increased eclipse depths 
of planet candidates orbiting low-metallicity stars. 

I now make a speculative proposition, which is that 
low-metallicity stars may host a higher proportion of 
planets that formed by gravitational instability than 
high-metallicity stars. While dust is the necessary raw 
ingredient for the protoplanetary cores that nucleate the 
growth of massive gaseous atmospheres in the core accre- 
tion scenario, dust is the enemy of planet formation by 
gravitational instability because it inc reases disk opac - 
ity, reducing the disk's ability to cool (|Cai et al.l f2006). 
Though planets that form by core accretion tend to be 
metal-rich even relative to their host stars an d by ne- 
cessity have solid cores (Miller fc Fortne"vll2011l ). planets 
that form by gravitational instability most likely have 
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the stella r composition and tend to b e less dense (see, 
however. iHelled fe~ Bodcnh eimerl [20Tol) . The idea that 
the giant planets orbiting low-metallicity stars predom- 
inantly formed by gravitational instability is interesting 
and warrants further thought. The theory that high- 
metallicity stars and low-metallicity stars form their gas 
giants b y different m e chanis ms is consistent with the re - 
sults of ISantos et~aT1 (|2004D and lUdrv fc Santos! (|2007l ). 
who found that the frequency of giant planets as a func- 
tion of [Fc/H] is flat at subsolar metallicities. If the 
dominant planet formation mechanism were not chang- 
ing with metallicity, one would expect a monotonically 
increasing planet detection rate with [Fe/H] for the entire 
metallicity range surveyed. 

While the eclipse depth-metallicity trend is not yet on 
firm statistical footing, it is an intriguing possibility that 
may provide insight into the Galactic history of planet 
formation. This work provides the first piece of evidence 
that giant planet structure, not just detection rate, de- 
pends on stellar metallicity. The planet populations of 



metal-rich stars and metal-poor stars may be very dif- 
ferent. It is even possible that differences in planetary 
structure with host star metallicity could indicate not 
only changes in the raw ingredients available for planet 
formation as the Galaxy becomes chemically enriched, 
but changes in the planet formation mechanism itself — 
planetary evolution across cosmic time. 
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Fig. 2. — Fraction of giant planet candidates with R p /R t > 0.13. Here the planets orbiting the lowest-metallicity hosts appear to be 
significantly bigger, relative to their parent stars, than the planets orbiting Solar-metallicity hosts. 




Fig. 3. — Left: A plot of [M/H] from the Kepler Input Catalog (KIC) as a function of spectroscopially determined [M/H] for published 
Kepler discoveries provides insight into the precision and accuracy of KIC metallicities. The best-fit line (black) is y = 0.88a; — 0.14, 
indicating that KIC metallicities are systematically lower than their spectroscopic counterparts. The precision of the KIC metallicities, 
CT[m/h] = 0-22, is indicated by the dash-dotted lines. The two outliers, Kepler-10 and Kepler-19, are both Solar-type, single dwarf stars. 
Right: Here I investigate whether a [M/H] is a function of T e ff in the selected sample, a relationship that could introduce statistical biases 
due to changes in both planet and companion star occurrence rates as a function of stellar mass. The solid line shows the mean [M/H] 
of the sample, while the dash-dotted line shows the best-fit line representing [M/H] vs. T c ff . The Spearman rank correlation coefficient is 
-0.19, indicating a weak dependence of metallicity on temperature. 
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Fig. 4. — Left: The Spearman rank correlation coefficient for planet radius as a function of log(g) is 0.04, with significance 0.55, indicating 
no trend that could influence sample selection. Right: [M/H] and log(g) values in the KIC are weakly correlated, with a Spearman rank 
coefficient of -0.13. However, my analysis indicates that KIC [M/H] is more strongly related to spectroscopic [Fe/H] than any other 
parameter. The dash-dot line shows the best fit representing [M/H] as a function of log(g). 




Fig. 5. — Monte Carlo simulations that take into account the errors in KIC metallicities suggest that there is a relationship between gas 
giant size and host star [M/H]. Here I show the fraction of systems with R v /R* > 0.13 as a function of host star [Fe/H] from 2 X 10 4 
Monte Carlo simulations, assuming errors on KIC metallicities are Gaussian and <T[m/h] = 0.22. The higher occurrence rate of large planets 
around the lowest-metallicity hosts, compared with Solar-metallicity and supersolar hosts, appears to be statistically significant. 
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TABLE 1 Giant planet candidates used in analysis 



KOI number 1 


KIC number 


Rp/R* 


Rp/R* error 


KIC [M/H] 


K00005.01 


8554498 


0.037 


0.0002 


0.12 


K00010.01 


6922244 


0.091 


0.0007 


0.00 


K00012.01 


5812701 


0.087 


0.0001 


-0.04 


K00020.01 


11804465 


0.117 


0.0000 


-0.16 


K00022.01 


9631995 


0.092 


0.0001 


0.29 


K00063.01 


11554435 


0.057 


0.0010 


0.12 


K00064.01 


7051180 


0.040 


0.0003 


-0.34 


K00082.01 


10187017 


0.034 


0.0025 


0.20 


K00089.02 


8056665 


0.022 


0.0006 


0.06 


K 00094 01 


6462863 

\J rt \J £t {J \J r.J 


0.070 


0.0024 


-0.76 


K000Q4 OS 


6462863 


0.038 


0.0016 


-0.76 


K00097.01 


5780885 


0.078 


0.0001 


0.05 


K00098.01 


10264660 


0.053 


0.0094 


-0.12 


K00100 01 


4055765 


0.045 


0005 


-0.37 


K00102.01 


8456679 


0.030 


0.0024 


-0.36 


K00105.01 


8711794 


0.039 


0.0006 


-1.08 


K00113.01 


2306756 


0.066 


0.0210 


0.10 


K00127.01 


8359498 


0.097 


0.0001 


0.17 


K00128.01 


11359879 


0.101 


0.0005 


0.36 


K00131.01 


7778437 


0.075 


0.0001 


0.16 


K00135.01 


9818381 


0.080 


0.0074 


0.18 


K00137.01 


8644288 


0.043 


0.0002 


0.17 


K001 37 02 


8644288 


0.062 


0.0013 


0.17 


K00138.01 


8506766 


0.094 


0.0006 


-0.22 


K001 39 01 


8559644 


0.058 


0.0022 


0.01 


K00144.01 


4180280 


0.035 


0.0047 


0.24 


K00161.01 


5084942 


0.029 


0068 


0.26 


K00183 01 


9651668 


0.122 


0.0190 


-0.14 


K00186 01 


12019440 


0.118 


0.0001 


0.02 


K00187 01 


7023960 


0.142 


0.0130 


0.10 


K00188.01 


5357901 


0.108 


0.0001 


0.26 


K00189.01 


11391018 


0.133 


0.0008 


-0.05 


K00190.01 


5771719 


0.112 


0.0006 


0.22 


K00191.01 


5972334 


0.115 


0.0011 


-0.19 


K00192.01 


7950644 


0.090 


0.0210 


-0.23 


K00193.01 


10799735 


0.129 


0.0002 


-0.12 


K00194.01 


10904857 


0.134 


0.0037 


-0.09 


K00195.01 


11502867 


0.114 


0.0019 


-0.19 


K00196.01 


9410930 


0.097 


0.0009 


0.10 


K00197 01 


2987027 


0.091 


0.0001 


0.00 


K00199.01 


10019708 


0.093 


0.0170 


0.10 


K00200.01 


6046540 


0.083 


0.0001 


0.15 


K00201.01 


6849046 


0.072 


0.0001 


0.19 


K00202.01 


7877496 


0.102 


0.0005 


0.12 


K00203.01 


10619192 


0.130 


0002 


0.39 


K00204.01 


9305831 


0.076 


0.0001 


-0.10 


K00205.01 


7046804 


0.092 


0.0020 


-0.17 


K00206.01 


5728139 


0.063 


0.0001 


-0.13 


K00208.01 


3762468 


0.089 


0.0007 


0.21 


K00209.01 


10723750 


0.069 


0.0002 


-0.04 


K00211.01 


10656508 


0.081 


0.0960 


0.01 


K00212.01 


6300348 


0.064 


0.0140 


-0.24 


K00214.01 


11046458 


0.111 


0.0035 


0.02 


K00216.01 


6152974 


0.065 


0.0004 


0.18 


K00217.01 


9595827 


0.134 


0.0001 


0.22 


K00229.01 


3847907 


0.049 


0.0002 


0.22 


K00242.01 


3642741 


0.056 


0002 


0.15 


K00254.01 


5794240 


0.184 


0.0012 


0.23 


K00256.01 


11548140 


0.123 


0023 


0.58 


K00261.01 


5383248 


0.027 


0.0067 


0.02 


K00345.01 


11074541 


0.035 


0.0063 


0.26 


K00348 01 


11194032 


0.038 


0003 


0.37 


K00351.01 


11442793 


0.083 


0.0064 


-0.11 


K00351.02 


11442793 


0.059 


0.0006 


-0.11 


K00353.01 


11566064 


0.064 


0.0012 


-0.09 


K00356.01 


11624249 


0.033 


0.0071 


-0.50 


K00367.01 


4815520 


0.044 


0.0085 


-1.53 


K00368.01 


6603043 


0.085 


0.0001 


-0.03 


K00372.01 


6471021 


0.081 


0.0093 


-0.56 


K00375.01 


12356617 


0.077 




-0.13 


K00377.01 


3323887 


0.078 


0.0030 


0.17 


K00377.02 


3323887 


0.084 


0.0017 


0.17 


K00398.01 


9946525 


0.092 


0.0036 


0.14 


K00401.01 


3217264 


0.041 


0.0003 


0.15 


K00401.02 


3217264 


0.043 


0.0034 


0.15 


K00410.01 


5449777 


0.102 


0.0050 


-0.40 


K00412.01 


5683743 


0.053 


0.0002 


-0.01 
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TABLE 1 Giant planet candidates used in analysis 



KOI number 1 


KIC number 


Rp/R* 


Rp/R* error 


KIC [M/H] 


K00415.01 


6289650 


0.062 


0.0630 


-0.33 


K00417.01 


6879865 


0.097 


0.0017 


-0.14 


K00418.01 


7975727 


0.115 


0.0004 


-0.32 


K00419.01 


8219673 


0.091 


0.0004 


-0.13 


K00421.01 


9115800 


0.115 


0.0002 


-0.07 


K00422.01 


9214713 


0.138 




-0.30 


K00423.01 


9478990 


0.085 


0.0003 


-0.18 


K00425.01 


9967884 


0.133 


0.0140 


0.07 


K00428.01 


10418224 


0.056 


0.0001 


-0.14 


K00433.01 


10937029 


0.049 


0.0180 


0.38 


K00433.02 


10937029 


0.113 


0.0013 


0.38 


K00458.01 


7504328 


0.077 


0.0062 


-0.20 


K00464.01 


8890783 


0.067 


0.0003 


0.17 


K00469 01 


9703198 


0.061 


0.0029 


-0.07 


K00523.01 


8806123 


0.063 


0.0015 


-0.09 


K00552.01 


5122112 


0.097 


0.0013 


-0.47 


K00554.01 


5443837 


0.069 


0.0034 


-0.08 


K00607.01 


5441980 


0.075 


0.0009 


-0.51 


K00609.01 


5608566 


0.089 


0.0110 


-0.08 


K00611.01 


6309763 


0.073 


0.0004 


-0.13 


K00617.01 


9846086 


0.177 


0.0210 


0.10 


K00620.01 


11773022 


0.072 


0.0009 


-0.08 


K00622.01 


12417486 


0.073 


0.0025 


-0.05 


K00625.01 


4449034 


0.062 


0.0097 


-0.06 


K00633 01 


4841374 


0.028 


0023 


-0.37 


K00674.01 


7277317 


0.038 


0.0042 


0.08 


K00680 01 


7529266 


0.060 


0.0001 


-0.46 


K 00684 01 


7730747 


0.041 


0.0042 


-0.01 


K00686.01 


7906882 


0.108 


0.0032 


-0.18 


K00716.01 


9846348 


0.061 


0.0027 


-0.33 


K00725.01 


10068383 


0.083 


0.0024 


-0.20 


K00728.01 


10221013 


0.099 


0.0010 


0.04 


K00737.01 


10345478 


0.064 


0.0020 


-0.28 


K00741.01 


10418797 


0.242 


0.0062 


0.06 


K00743.01 


10464078 


0.087 


0.0210 


-0.56 


K00745.01 


10485250 


0.092 


0.0006 


-0.20 


K00753.01 


10811496 


0.102 


0.0072 


-0.13 


K00760.01 


11138155 


0.107 


0.0008 


0.01 


K00763.01 


11242721 


0.110 


0.0009 


-0.13 


K00764.01 


11304958 


0.047 


0.0004 


0.23 


K00767.01 


11414511 


0.128 


0.0008 


0.02 


K00771.01 


11465813 


0.124 




-0.11 


K00772.01 


11493732 


0.070 


0.0044 


0.00 


K00774.01 


11656840 


0.143 


0.0004 


-0.31 


K00779 01 


11909839 

_1_ JU^C 


0.109 


0.0210 


-0.55 


K00782.01 


11960862 


0.048 


0.0930 


-0.31 


K00791 01 


12644822 


0.071 


0003 


-0.01 


K00797.01 


3115833 


0.077 


0.0007 


-0.52 


K00801.01 


3351888 


0.081 


0.0130 


0.18 


K00802.01 


3453214 


0.135 


0.0049 


-0.41 


K00805.01 


3734868 


0.119 


0.0095 


-0.10 


K00806.01 


3832474 


0.093 


0.0073 


-0.10 


K00806.02 


3832474 


0.125 


0.0005 


-0.10 


K00809.01 


3935914 


0.114 


0.0650 


-0.39 


K00813.01 


4275191 


0.085 


0.0003 


-0.28 


K00815.01 


4544670 


0.101 


0.0130 


-0.16 


K00822.01 


5077629 


0.128 


0.0015 


-0.40 


K00823.01 


5115978 


0.075 


0.0015 


-0.43 


K00824.01 


5164255 


0.122 


0.0013 


-0.21 


K00830.01 


5358624 


0.134 


0.0001 


0.15 


K00838.01 


5534814 


0.072 


0.0010 


-0.10 


K00840.01 


5651104 


0.096 


0.0031 


-0.09 


K00843.01 


5881688 


0.053 


0.0049 


0.20 


K00846.01 


6061119 


0.166 


0.0008 


-0.52 


K00847 01 


6191521 


0.054 


0330 


-0.57 


K00850.01 


6291653 


0.092 


0.0042 


-0.19 


K00851.01 


6392727 


0.056 


0.0007 


-0.12 


K00855.01 


6522242 


0.138 


0.0002 


-0.46 


K00856.01 


6526710 


0.140 


0.0025 


-0.15 


K00858.01 


6599919 


0.091 


0.0027 


-0.25 


K00865.01 


6862328 


0.074 


0.0360 


-0.09 


K00868.01 


6867155 


0.161 


0.0025 


-0.71 


K00871.01 


7031517 


0.208 


0.0044 


-0.40 


K00872.01 


7109675 


0.084 


0.0027 


-0.11 


K00876.01 


7270230 


0.144 


0.0066 


-0.12 


K00878.01 


7303253 


0.041 


0.0031 


-0.27 


K00880.02 


7366258 


0.055 


0.0004 


-0.01 
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Dodson- Robinson 
TABLE 1 Giant planet candidates used in analysis 



KOI number 1 


KIC number 


Rp/R* 


Rp/R* error 


KIC [M/H] 


K00882.01 


7377033 


0.151 


0.0015 


-1.24 


K00883.01 


7380537 


0.167 


0.0005 


-0.37 


K00889.01 


757450 


0.114 


0.0003 


-0.09 


K00890.01 


7585481 


0.077 


0.0003 


0.03 


K00895.01 


7767559 


0.107 


0.0011 


-0.26 


K00897.01 


7849854 


0.109 


0.0001 


0.27 


K00902.01 


8018547 


0.080 


0.0007 


-0.50 


K00903 (11 


8039892 


0.077 


0002 

v / . \ / I/O- 


-0.32 


K00908.01 


8255887 


0.081 


0.0002 


0.13 


K00913.01 


8544996 


0.122 


0.0049 


-0.28 


K00918 01 


8672910 


0.111 


0003 


-0.07 


K00929 01 


9141746 


0.078 


0002 


0.14 


K00931.01 


9166862 


0.116 


0.0075 


0.32 


K00941.01 


9480189 


0.043 


0005 


0.32 


K00941 03 


9480189 


0.052 


0.0120 


0.32 


K00951.01 


9775938 


0.046 


0.0120 


0.41 


K00956.01 


9875711 


0.044 


0.0170 


0.39 


K00960.01 


8176650 


0.183 


0.0005 


-0.09 


K00972.01 


11013201 


0.019 


0.0025 


-0.06 


K00981.01 


8607720 


0.013 


0.0012 


-0.11 


K00988.01 


2302548 


0.033 


0.0021 


0.14 


K01003.01 


2438502 


0.141 


0.0007 


-0.09 


K01005.01 


5780460 


0.062 


0.0430 


-0.21 


K01089.01 


3247268 


0.083 


0.0003 


-0.13 


K01089.02 


3247268 


0.049 


0.0060 


-0.13 


K01159.01 


10354039 


0.054 


0.0039 


0.29 


KO1 1 76 01 


3749365 


0.157 


0003 


0.38 


K01177.01 


3547091 


0.130 


0.0030 


-0.57 


K01193.01 


3942446 


0.106 


0.0380 


-0.53 


K01208.01 


3962440 


0.060 




-0.20 


K01221.02 


3640905 


0.012 


0.0059 


0.42 


K01227.01 


6629332 


0.120 


0.0013 


-0.20 


K01241.01 


6448890 


0.020 


0039 


0.16 


K01241.02 


6448890 


0.013 


0.0051 


0.16 


K01242.01 


6607447 


0.058 


0006 


0.00 


K01257.01 


8751933 


0.080 


0.0005 


-0.18 


K01261.01 


8678594 


0.064 


0.0005 


0.12 


K01268.01 


881 3698 


0.074 




-0.07 


K01285.01 


10599397 


0.081 


0.0029 


-0.09 


K01288.01 


10790387 


0.084 


0.0005 


-0.20 


K01299.01 


10864656 


0.026 


0.0002 


-0.24 


K01335.01 


4155328 


0.040 


0.0210 


-0.28 


K01353.01 


7303287 


0.102 


0.0011 


-0.08 


K01385.01 


9278553 


0.198 


0.0003 


-0.32 


K01391.01 


8958035 


0.079 


0.0016 


-0.21 


K01419.01 


11125936 


0.053 


0.0035 


-0.29 


K01426.02 


11122894 


0.067 


0.0028 


-0.12 


K01459.01 


9761199 


0.075 


0.0010 


0.10 


K01474.01 


12365184 


0.062 


0.0005 


-0.33 


K01477.01 


7811397 


0.122 




0.00 


K01486.01 


7898352 


0.094 




-0.24 


K01540.01 


5649956 


0.197 


0.0011 


-0.44 


K01543.01 


5270698 


0.148 


0.0003 


-0.24 


K01546.01 


5475431 


0.101 


0.0008 


0.12 


K01553.01 


7951018 


0.069 


0.0006 


0.21 


K01557.01 


5371776 


0.039 


0.0350 


0.52 


K01561.01 


4940438 


0.052 


0.0016 


-0.60 


K01574.01 


10028792 


0.062 


0.0003 


0.05 


K01587.01 


9932970 


0.080 


0.0067 


-0.08 



1 Kepler Object of Interest number 



Eclipse depth-metallicity correlation 
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TABLE 2 

Kepler DISCOVERIES USED FOR METALLICITY ACCURACY ANALYSIS 



Name 



KIC number Spectroscopic [Fe/H] 



Reference 





1 1 44R443 


—0.15 ± 0.1 


\\ O T~\ 1 A. 


1 1 K^QO^ 


u.i ( in u.uu 


R'onlpr 7 
J. vl ]Jlv_ 1 i 


O / 0\J000 




IVCU1C1 o 


6999944 


_n nfi + n 03 

U.UU _l_ U.UU 






n i o _l n 04 

U. IC U.U'i 


K"pnlpr 10 

J. vl ]Jlv_ 1 ± \l 


1 1 Q041 M 


_n i c _l o r\A 

— U.-LO IC U.U'I 


T< O 1 £i 1 "1 


fi^zLI Q90 


U.UU IC U. J.U 


K"pnlpr 1 9 


1 1 80446^ 


n 07 4- 04 


TCor^lpr 1 A 


1 09fi4fi(S0 


n i o _|_ n r\f. 


K OT^lor 1 ^ 




n Qfi -U 07 
u.ou zc u.u/ 


TCpnlpr IfiAR 


1 9fi447fiQ 


—0.3 i 0.2 


R'pnlpr 1 7 


1 Ofil Q1 Q9 


u. ^u m u. J. 


K"pnlpr 1 8 


£644988 


n 9 4- 04 

U. ^ IC U.U'i 


Kpnlpr-1 Q 


2571 238 


U.-LU _l_ U.UU 




uoouou^t 


09 4- 04 
u.u^ zc u.u^t 


Kepler-22 


10593626 


-0.29 ± 0.06 


Kepler-23 


11512246 


-0.09 ±0.14 


Kepler-26 


9757613 


-0.21 ±0.08 


Kepler-27 


5792202 


0.41 ± 0.04 


Kepler-30 


3832474 


0.18 ±0.027 


Kepler-34AB 


8572936 


-0.07 ±0.15 


Kepler-35AB 


9837578 


-0.34 ±0.2 


KOI-135 


9818381 


0.33 ±0.11 


KOI-196 


9410930 


-0.09 ±0.16 


KOI-204 


9305831 


0.26 ±0.1 


KOI-254 


5794249 


0.13 ±0.13 


KOI-423 


9478990 


-0.29 ±0.1 




6185331 


0.11 ±0.15 



Torres ct al. (2008) 
Borucki et al. (2010) 
Latham et al. (2010) 



Jenkins et al. (2010) 
Holman et al. (2010) 
Batalha et al. (2011) 
Lissauer^e^il^ (2011) 
Fortnev et al. (2011) 
Buchhave et al. (5uTl) 
Endl et al~720lT) 

Doyle et al. 12011) 
Bonomo et al. (2011) 
Cochran et~aL (2011) 

Ballard et al. (2011) 

Fressin et al. (2011) 
Borucki et al. (_2012b) 
Ford et al. (2012) 

Steffen et al.(2uT2) 

Steffen et a~ (20r2) 
Fabrvckv et al (2012) 

Welsh et1T720i*2) 

Welsh et al. (2U12*) 
Bonomo et al. (2011) 
SanterrreetaL 
Bonomoetal. 
JohnsonetaL 

Bouch^et_aL 

Fischer et al. 



!011) 
f207l) 

'mm 

' mm 

(2oTT) 
(2012) 



